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Abstract. - When two chemically passivated solids are brought into contact, interfacial 
interactions between the solids compete with intrabulk elastic forces. The relative importance of 
these interactions, which are length-scale dependent, will be estimated using scaling arguments. 
If elastic interactions dominate on all length scales, solids will move as essentially rigid objects. 
This would imply superlow kinetic friction in UHV, provided wear was absent. The results of 
the scaling study depend on the symmetry of the surfaces and the dimensionalities of interface 
and solids. Some examples are discussed explicitly such as contacts between disordered three- 
dimensional solids and linear bearings realized from multiwall carbon nanotubes. 



Introduction. - Many small-scale devices cannot be miniaturized further, because friction 
and wear appear to be exceedingly large in nanoscale machines. This limitation seems to arise, 
because the surface to volume ratio is large in small systems and thus surface forces such as 
friction become relatively large. However, many theoretical studies and most 

recently an increasing number of experiments indicate that shear forces can 

be tremendously small between two atomically flat surfaces. These findings spur the hope for 
new avenues to reduce friction in nano-scale applications. 

Superlow friction was first suggested by Hirano and Shinjo Their calculations in- 

dicated that no instabilities occur when two atomically smooth copper solids slide past each 
other provided that they are sufficiently misaligned. The absence of instabilities in such con- 
tacts implies that kinetic friction approaches zero at infinitely small velocities, even when 
thermal fluctuations are absent. Hirano and Shinjo called such superlow friction superlubric- 
ity. This term is often considered unfortunate, because one may expect zero friction in the 
sliding state in analogy to superconductivity and superfluidity. However, due to the emission 
of sound waves, which also occur when instabilities are absent, a drag force linear in velocity 
will persist. It might therefore be more appropriate to call the effect structural lubricity, as 
the low friction arises from the structural incompatibility of the two contacting solids rather 
than from a Bose Einstein condensate of bosons. The main reason for the low friction is that 
lateral forces between two non-matching, rigid solids cancel systematically for incommensu- 
rate interfaces and stochastically between disordered surfaces, so that the average force per 
unit area decreases quickly with increasing contact area. 
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In principle, two solids can pin even when their intrinsic surface corrugations do not match. 
As suggested already by Amontons in 1699 |12|. solids can deform elastically so that the de- 
formed structures eventually interlock geometrically. However, elastic deformations would 
only produce the typically observed friction laws, if Prandtl's condition for (elastic) pinning 
was satisfied, namely the condition that for a given macroscopic configuration (i.e. a given rel- 
ative displacement of the two solids), several non-equivalent, mechanically stable, microscopic 
configurations exist ^Sj- In that case, externally imposed sliding would lead to instabilities 
(or plucking motion) and ultimately to dissipated energy. This mechanism is discussed in 
more detail elsewhere ^31- I n the Prandtl (or Tomlinson) model, multistability and conse- 
quently clastic hysteresis and finite kinitec friction occur when the maximum curvature of 
the substrate potential exceeds the spring constant with which the surface atom is bound to 
its lattice site. Although the Prandtl Tomlinson model is only one dimensional, it shows the 
importance of the competition between intrabulk clastic restoring forces and interfacial inter- 
actions. While there is a large literature on such competitions in the field of vortex motion in 
superconductors |15| . no generic, simple, and quantitative analysis is known to the author, in 
which the effect of interfacial order and dimensionality d- m t as well as the solids' dimension- 
ality d Q bj is discussed. In the context of mechanical friction between solids, approaches are 
either computational, highly complex, and/or the issue of dimensionality and orientational 
alignment is discussed qualitatively. 

The subject of this letter is a generic, quantitative analysis of how fast interfacial inter- 
actions and intrabulk clastic forces increase as the length scale (either of the object or of the 
description) increases. In this analysis, we will assume that it is legitimate to treat a block 
of length C as effectively rigid as long as the intrabulk stiffness (to be defined below) on that 
length scale is greater than the interfacial stiffness on the same length scale. The reason is that 
Prandtl's multi-stability criterion cannot be satisfied in that situation. We will then analyze 
whether this assumption still holds when the length scale is increased ever further. If, at some 
point, the rigid-block assumption breaks down, then some fractions of the block may be able 
to move in a partially uncorrelated fashion, i.e., multistability and consequently instabilities 
can occur that do not incorporate the motion of the whole block. In the latter case, the 
interface would not show structural lubricity. It has to be emphasized that the present study 
cannot predict with certainty whether an interface really loses its structural lubricity when 
the rigid block assumption breaks down. Our discussion would only suggest that there is the 
possibility for the loss of structural lubricity, or - in Aubry's words |16II17| - for the possibility 
of the breaking of analyticity. In the remaining part of this letter, I will present the scaling 
arguments for interfacial and intrabulk interactions. These arguments will then be applied to 
selected cases. 

Theory. - The starting point of our analysis is displayed in fig. 1(a). Two solids with flat 
interfaces are brought into contact. Let us assume that the interatomic interactions within 
each solid are strong (covalent, metallic, ionic bonds or a combination thereof), while the 
interactions between the solids are physical, i.e., relatively weak Lennard Jones- type forces 
between chemically passivated solids. In such a situation, one may assume that one can 
approximate the intrabulk interactions with harmonic springs of stiffness &b°uik> wnnc (except 
for small corrections due to small atomic displacements) the interactions of surface atoms with 
the opposing surface (called the substrate) can be described by a function that has the same 
periodicity a s as the substrate, i.e., in one dimension V(x) = Vo cos(27rx/a s ), see Fig. 1(b). 
The maximum curvature of the substrate potential will be denoted as atomic-scale interfacial 
stiffness k\°%. It is possible to estimate the order of magnitude for the values of atomic-scale 
intrabulk and interfacial stiffness, k]^ ik and k\°^, respectively. Their dimensions are given by 
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Fig. 1 - Schematic view of the interaction between chemically passivated solids. Part (a) shows a 
snapshot of a simulation (whose results have been reported in Ref. 5 1) between two solids. Part (b) 
represents the coupling of a surface atom to its neighbours (reflected by the harmonic springs) and to 
the substrate (reflected by the sinusoidal line). The parabola indicates the maximum curvature of the 
atom-substrate potential, which is called the atomic-scale interfacial stiffness fc}"'. Part (c) describes 
the scaling procedure for a ID elastic chain. 



elastic constant times interatomic spacing a. Strong bonds lead to bulk modules of typically 
more than 40 GPa, while physical bonds are usually responsible for bulk modules in the 
order of 4 GPa. a can be roughly estimated with 2 A for chemical bonds and with 3 A for 
physical bonds. Based on these orders of magnitude, one may argue that fc^° u c lk w 40 GPa 2 A 
and k\°% w 4 GPa 3 A. It is therefore obvious in such a situation that on the atomic scale, 
each (surface) atom has exactly one well-defined mechanical equilibrium position, once the 
coordinates of its intrabulk neighbours relative to the substrate are known. The situation is 
similar to that in the Prandtl or Tomlinson model when the elastic coupling to the lattice 
site is greater than the maximum curvature of the substrate potential. Thus, atomic-scale 
instabilities cannot occur, which, however, does not preclude instabilities on larger length 
scales. 

In the following, we will increase the length scale of our description from C = a to C = 2a. 
Since fcbuik i s much greater than font on C = a, let us group together linear blocks of length 
2cr and treat them as rigid objects on the new length scale. Of course, the values for fcbuik 
and fci n t have to be renormalized for the new length scale. In the following renormalization, 
we will treat intrabulk and interfacial interactions completely separately. This is certainly 
a simplistic procedure as it suppresses generalized terms in the Hamiltonian that arise in a 
true renormalization group (RG) treatment of special case studies ^8]. Thus, the present 
examination does not attempt to be a rigorous RG study. However, it contains the leading 
contributions in the first few renormalization steps and, more importantly, it illustrates in 
an easy-to-digest fashion relevant key features for more detailed studies of specific systems. 
Moreover, the trends that will be obtained at the end of this analysis are in very good agree- 
ment with many computational and theoretical studies. The reason for this agreement may be 
that the neglected additional terms are irrelevant as long as the rigid-body hypothesis remains 
valid in the present treatment. 
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The coarse-graining procedure for the elastic interaction is shown in fig. 1(c) for a ID 
elastic solid. Each time we group together two atoms into one new 'coarse-grained atom' by 
doubling the lattice constant in our description, we have to decrease the elastic stiffness of the 
springs by a factor of two, in order to maintain the same continuum limit of the chain's elastic 
constants. Thus, the effective elastic stiffness of an elastic chain on a length scale C (or simply 
the stiffness of a linear chain of length C) is given by fct>uik(£) = ^buik/^- ^ tne elastic object 
is two dimensional, <i bj = 2, springs are coupled in series and in parallel. Therefore, fcbuik(£) 
remains unaltered when the scale is changed in 2D. Each additional dimension increases the 
effect of serial coupling and therefore 

fcbuikOC) = (W^Xuik- (1) 

As a short side comment, we note that the dimension d Q bj = 2 (at which solids neither 
become harder nor softer as their length scales increase) is identical to the dimension above 
which the Debye Waller factor of a (topologically) perfect crystal remains finite, which implies 
that Bragg peaks can be observed for c? bj > 2 even when thermal fluctuations are present, 
while they are unstable against thermal fluctuations for d bj < 2 

The coarse-graining procedure for the interfacial interactions depends on the symmetry of 
the surfaces. Let us first consider two commensurate solids. Each surface atom in the slider 
experiences exactly the same curvature in the substrate potential as every other surface atom. 
Therefore, if we group together C dobi atoms into one coarse-grained atom, C dint surface atoms 
will give a coherent contribution to the interfacial stiffness, leading to 

= (£/(7)*»'Afe (2) 

where fc?° t m (£) is the interfacial stiffness of commensurate solids on a length scale C. 

If at least one of the two surfaces is disordered, then the substrate potential that individual 
surface atoms experience will be random. The average surface curvature of interfacial stiffness 
has to be zero for flat surfaces and the variance be well defined. For the sake of simplicity, 
but without loss of generality, we will consider small correlation lengths. (Large correlation 
lengths lead to larger pre-factors in the following equation, but leave the scaling unaltered.) 
If we sum up the individual (random) contributions fc; n t, then the coarse-grained interfacial 
stiffness kf^(C) will also be random with zero mean and a variance that scales according to 
the laws of large numbers as 0|: 

({k^(C)} 2 )^(C/aY^((k^) 2 ). (3) 

If both surfaces are incommensurate, then the annihilation of lateral surface forces that 
individual surface atoms experience is rather systematic. As a consequence, we find that 
(L) does not increase systematically with £, hence 

ki»m « (4) 

This equation might be surprising, however, it is easily verified exemplarily for two simple 
cubic lattices of identical alignment but different lattice constants. It is important to mention 
that the (suppressed) proportionality factor in the previous two equations will in general 
depend on the length scale C. However, these prefactors will not increase systematically with 
L and for each specific case, it will be possible to define an upper bound for the proportionality 
coefficient. This upper bound will mainly be based on geometric considerations such as relative 
orientation of the solids and lattice misfit. For instance, the upper bound of the proportionality 
coefficient in eq. (@J will become large close to commensurability. 
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As discussed in the introduction, the crucial step is to compare how quickly fc; nt and fctmik 
change with the (coarse-grain) length C. If the reduced interfacial stiffness k{L) defined as 
the ratio of interfacial and intrabulk stiffness 

n{C) = k int (£)/k hulk {£) (5) 

remains smaller than unity for every value of C, then the rigid body assumption remains valid. 
The reason is that there is a well-defined equilibrium position of a coarse-grained atom for 
every given coordinates of neighbours on every length scale. Consequently, the system shows 
structural lubricity according to Prandtl's condition. If, in the scaling limit C — > oo or for any 
intermediate value of C, the rigid body hypothesis potentially breaks down, i.e., k(C) > 1, 
then the system might lose its structural lubricity. To summarize, our analysis only addresses 
the question whether or not we have to expect the solids to move as rigid units. When they 
do not move as a rigid unit, our present treatment cannot state whether motion can occur in 
terms of (continuous) solitons, in which case even athermal kinetic friction would vanish at 
zero sliding velocity, or in terms of (discrete) dislocations with finite kinetic friction. 

Examples. One-dimensional, incommensurate solids. The Frenkel-Kontorova (FK) 
model, which is a ID elastic chain embedded in a sinusoidal external potential, is the generic 
model for this class |2QI21j . The FK model is one of the most studied models in the context of 
friction. We can easily see from eqs. |0J, and JSJ that k increases linearly with C in this 
class. Thus, from our discussion, we have to allow for the possibility that the system shows 
friction. However, as emphasized before, a more detailed study has to be done whenever k 
exceeds unity. Systematically correct solutions of the FK model were provided by Aubry and 
coworkers, in particular Peyrard |16U17j . They found that for every irrational ratio a = a c /a s 
of the chain's average lattice constant a c and the substrate's lattice constant a s , there is a 
critical ratio of n c = fcmt/fcbuik above which analyticity is broken, or structural lubricity is 
lost. k c is close to unity when a equals the golden mean and smaller for other values of a. 

Nanotubes are also quasi ID structures that can form ID-incommensurate interfaces, for 
instance when a nanotube is placed onto graphite out of registry or when two nanotubes are 
nested into each other in an incommensurate fashion (i.e. in the armchair/zigzag configura- 
tion). Thus nanotubes do not necessarily behave structurally lubric, even when out of registry. 
We yet have to expect superlow friction, even structural lubricity, because the reduced inter- 
facial stiffness n at the atomic scale will be much smaller than unity. An estimate would 
indicate that fcbuik is proportional to the stiffness of graphite's in-plane covalent bond times 
the diameter R of the nanotube, while fc; n t would be a typical value for fcj nt based on Lenard 
Jones type interactions, hence k(ij) w C44CT jC\\R < 10 -3 , where C\\ and C44 are elastic 
constants of graphite. For nested nanotubes, one may expect even smaller values of Kmt(c) 
due to periodic-boundary conditions of the toruses. Thus, while nanotubes will probably not 
move as rigid blocks, the value of n(a) is so small that analyticity will most likely not be 
broken. 

Of course, this estimate docs not include effects due to the opening of the nanotubes. 
Recent computer simulations indicate that dangling bonds are present at the entrance of 
multiwall nanotubes. These bonds are chemically active and explain discrepancies between 
theoretically predicted and experimentally measured shear forces [2~2"] . 

Three-dimensional disordered solids with two-dimensional interface. A generic case of solid 
sliding is that between two 3D solids with a 2D interface. If the surface is disordered, then 
the interfacial stiffness increases as quickly with length scale as the intrabulk elastic forces. 
In situations, where restoring forces show the same scaling as random forces, the fluctuations 
typically win on exponentially large length scales. (This phenomenon is well-known when 
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studying elastic neutron scattering of 2D solids at non-zero temperature, where the intensities 
of Bragg reflexes decrease logarithmically with increasing system size [E].) Indeed, Pers- 
son [23] and Sokollof [25 found such logarithmic laws. They studied the elastic coherence 
length A (also called Larkin length |T2], the length scale below which a block is believed to 
move as a highly correlated unit and above which it moves less correlated) , and found that A 
depends exponentially on elastic constants and other parameters. Thus, the present scaling 
study is in qualitative agreement with their results for sufficiently flat, disordered interfaces. 
If, however, the surfaces are fractal, which is the case for almost any macroscopic object, then 
the random lateral forces will cancel less quickly with increasing coarse-graining than for flat 
surfaces. It may therefore be more realistic to expect a powerlaw dependence of A on elastic 
constants and atomic-scale interfacial stiffnesses. Yet, the net friction force that comes from 
pinning due to surface corrugation would remain extremely small. 

Three-dimensional, incommensurate crystals. The ratio of interfacial and intrabulk stiff- 
nesses becomes very small for 3D crystals that form an atomically smooth, incommensurate 
interface, i.e., k decreases with 1/yL. Therefore, if plastic flow and surface contaminations 
are absent, we should expect structural lubricity, unless the orientation of the like solids is 
close to commensurability. Clean, incommensurate interfaces have been studied exceedingly 
by computer simulation and essentially all studies based on interatomic potentials point to 
structural lubricity, in some cases even for non-passivated surfaces, i.e., the original studies 
on superlubricity 0[2] . A more detailed overview is given in Ref. ^5 • 

Conclusions and Outlook. - In this letter, I discussed the competition between interfacial 
interactions and intrabulk elastic forces when two chemically passivated solids are brought 
into mechanical contact. Simple scaling arguments were provided that allow one to estimate 
whether structural lubricity (essentially zero kinetic friction in the small-velocity limit) should 
be expected in UHV (as long as plastic deformation does not occur). Some specific examples 
were discussed in terms of the scaling arguments. For instance, the analysis showed that 
embedded nanotubes are not structurally lubric in a strict sense, even though interfacial 
interactions are too weak to break analyticity. Two incommensurate 3D crystals forming a 
flat 2D interface should almost always be structurally lubric. However, if the surfaces are 
disordered and moreover fractal, then structural lubricity will most likely be lost. The latter 
result would imply that different asperities in a typical multi-asperity contact do not move in 
a correlated way. 

The present manuscript focused on flat or at least self-similar surfaces. A separate treat- 
ment would have to be done for curved surfaces. First steps in that direction have been 
taken in a computer simulation study of an AFM tip that is pressed against a substrate |25| . 
Another interesting issue would be to study the effect of thermal and quantum fluctuations 
on the drag forces between structurally lubric solids. (Most studies so far have been done 
for the ID FK model, which probably does not show the same universal behaviour as 3D 
solids, as the reduced interfacial stiffness tends to infinity at large length scales for the ID 
FK model.) Studying the (quantum) friction between two incommensurate He 4 solids is cer- 
tainly computationally feasible. It may be speculated that true superlubricity occurs below 
the temperature at which solid He condenses into a 'supercrystal'. 

It will certainly remain an experimental challenge to prepare sufficiently flat, clean, and 
chemically inert surfaces to make technological use of structural lubricity. Moreover, structural 
lubricity will break down at normal pressures that are high enough to deform solids plastically. 
At least this latter point might not be an issue in nano-scale devices, as they are less rough 
than normal macroscopic surfaces and therefore have less extreme pressure distributions. 
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